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Abstract 

on : We have investigated an LRS Bianchi Type I models with bulk vis- 

cosity in the cosmological theory based on Lyra's geometry. A new class 
of exact solutions have been obtained by considering a time-dependent 
displacement field for a constant value of the deceleration parameter and 
viscosity coefficient of bulk viscous fluid is assumed to be a power function 
^ | of mass density. The physical behaviour of the models is also discussed. 
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1 Introduction 



After Einstein, in 1917, developed his general theory of relativity in which grav- 
itation is described in terms of geometry, Weyl, in 1918, proposed a more general 
i theory in which electromagnetism is also described geometrically. However, this 

theory, based on non-integrability of length transfer, had some unsatisfactory 
features and did not gain general acceptance. Later Lyra 1 suggested a modifi- 
cation of Riemannian geometry, which may also be considered as a modification 
of Weyl's geometry, by introducing a gauge function into the structureless man- 
ifold which removes the non-integrability condition of the length of a vector 
under parallel transport and a cosmological constant is naturally introduced 
from the geometry. In subsequent investigations, Sen |2j & Sen and Dunn [3] 
proposed a new scalar-tensor theory of gravitation and constructed an analog 
of the Einstein field equations based on Lyra's geometry. 

Halford^] pointed out that the constant displacement vector field </>,; in 
Lyra's geometry plays the role of a cosmological constant in the normal gen- 
eral relativistic treatment. Halford[S] showed that the scalar-tensor treatment 
based on Lyra's geometry predicts the same effects, within observational lim- 
its, as in Einstein's theory. Several authors (Bhamra[S], Karade and Borikar[7], 
Kalyanshetti and Wagmode|H], Reddy and InnaiahS], Beesham|10|. Reddy and 
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Venkateswarlu [TT] and Soleng^^l) have studied cosmological models based on 
Lyra's geometry with a constant displacement field vector. However, this re- 
striction of the displacement field to be a constant is a coincidence and there is 
no a priori reason for it. Singh and Singh[T3| — 115| and Singh and Desikan[TC| 
have studied Bianchi Type I, III, Kantowski-Sachs and a new class of models 
with a time dependent displacement field and have made a comparative study 
of Robertson- Walker models with a constant deceleration parameter in Ein- 
stein's theory with a cosmological term and in the cosmological theory based on 
Lyra's geometry. Recently Pradhan and Vishwakarma[T7| investigated a new 
class of an LRS Bianchi Type-I cosmological models in Lyra geometry. Though 
the displacement vector has no clear and unambiguous interpretations some ef- 
forts have been made to treat the constant displacement vector as the analogue 
of the cosmological constant. Soleng^^l nas pointed out that the cosmologies 
based on Lyra's manifold with a constant gauge vector <fr will either include 
a creation field and be equal to Hoyle's creation field cosmology [T^j — |2D] or 
contain a special vacuum field which together with the gauge vector may be 
considered as a cosmological term. In the latter case, the solutions are same 
as those of general relativistic cosmologies with a cosmological term. Recently 
Behnke et a2.[2D also pointed out to an alternative description of the new cos- 
mological supernova data without a A-term as evidence for Weyle's geometry 
of similarity, where Einstein's theory takes the form of the conformal-invariant 
theory of a massless scalar field|22| — 125). As it has been shown by Weyl al- 
ready in 1918, conformal-invariant theories correspond to the relative standard 
of measurement of a conformal-invariant ratio of two intervals, given in the ge- 
ometry of similarity as a manifold of Ricmannian geometries connected by 
conforms! transformations. This ratio depends on nine components of metrics 
whereas the tenth component became the scalar dilation field that cannot be 
removed by the choice of the gauge. In the current literature j57j this peculiarity 
of the conformal-invariant version of Einstein's dynamics has been overlooked. 
The energy constraint converts this dilation into a time-like classical evolution 
parameter which scales all masses including the Planck mass. In the confor- 
mal cosmology (CC), the evolution of the value of the massless dilation field 
(in the homogeneous approximation) corresponds to that of the scale factor in 
standard cosmology (SC). Thus, the CC is a field version of the Hoyle-Narlikar 
cosmology |28). where the redshift reflects the change of the atomic energy levels 
in the evolution process of the elementary particle masses determined by the 
scalar dilation field [23 EHHZ3- The CC describes the evolution of the confor- 
mal time, which has a dynamics different from that of the standard Friedmann 
model. Behnke et a/. [22 have also discussed as an observational argument 
in favour of the CC scenario that the Hubble diagram (effective magnitude- 
rcdshift-rclation: m(z)) including the recent SCP dataj^Hl can be described 
without a cosmological constant. 

Most studies in cosmology involve a perfect fluid. However, observed physi- 
cal phenomena such as the large entropy per baryon and the remarkable degree 
of isotropy of the cosmic microwave background radiation suggest that we should 
analyse dissipativc effects in cosmology. Furthermore, there are several processes 
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which are expected to give rise to viscous effects. These are the decoupling of 
neutrinos during the radiation era and the decoupling of radiation and matter 
during the recombination era. Bulk viscosity is associated with the GUT phase 
transition and string creation. The model studied by Murphv|38| possessed an 
interesting feature in that the big bang type of singularity of infinite space-time 
curvature does not occur a finite past. However, the relationship assumed by 
Murphy between the viscosity coefficient and the matter density is not accept- 
able at large density. The effect of bulk viscosity on the cosmological evolution 
has been investigated by a number of authors in the framework of general theory 
of relativitv|3U j — |4*5 | . This motivates us to study a cosmological viscous fluid 
model. 

The Einstein's field equations are a coupled system of highly nonlinear dif- 
ferential equations and we seek physical solutions to the field equations for their 
applications in cosmology and astrophysics. In order to solve the field equations 
we normally assume a form for the matter content or that spacetime admits 
killing vector symmetries 46 . Solutions to the field equations may also be gen- 
erated by applying a law of variation for Hubble's parameter which was proposed 
by BermanJI?]. It is interesting to obverse that the law yields a constant value 
for deceleration parameter (DP). The variation of Hubble's law as assumed is 
not inconsistent with observation and has the advantage of providing simple 
functional forms of the scale factor. In simplest case the Hubble law yields a 
constant value for the DP. The cosmological models with constant deceleration 
parameter (CDP) may be divided into two categories. The first category of 
models with CDP is that of models where the cosmic expansion is driven by big 
bang impulse; all the matter and radiation energy is proposed at the big bang 
epoch and the universe has started with singular origin. In the second cate- 
gory of models with CDP, the universe has a non-singular origin and the cosmic 
expansion is driven by the creation of matter particles. It is worth observing 
that most of the well-known models of Einstein's theory and Brans-Deke the- 
ory with curvature parameter k — 0, including inflationary models, are models 
with CDP. It also measures the deviation from linearity of growth of the scale 
factor. In earlier literature cosmological models with an CDP have been stud- 
ied by Berman|4*7|. Berman and Gomide0H|> Johri and Desikan fUfl I5U| . Singh 
and Dcsikan[TB|. Maharaj and Naidoo[5T]. Pradhan et aL |17lf4"Ull42| and others. 
This has provided us the motivation to study models with CDP. 

In this paper, we have investigated bulk viscous Locally Rotationally Sym- 
metric (LRS) Bianchi Type I cosmological models based on Lyra's geometry 
with a time dependent displacement field. It is remarkable to note here that the 
time dependent displacement vector may lead to the singularity free model 52 . 
We have obtained exact solutions of the field equations by assuming the decel- 
eration parameter to be constant. The physical behaviour of these models have 
been discussed. 
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2 Field Equations 



The metric for LRS Bianchi Type I spacetime is 

ds 2 = dt 2 - A 2 dx 2 - B 2 (dy 2 + dz 2 ), (1) 

where, A and B are functions of x and t. The energy momentum tensor in the 
presence of bulk stress has the form 

Tij = (P + P)uiUj ~ pg lj , (2) 

together with comoving coordinates u l Ui = 1 where Ui = (0, 0, 0, 1) and 

p = p + S,u l i and. (3) 

Here p,p,p,£, and u are, respectively, the energy density, isotropic pressure, 
effective pressure, bulk viscous coefficient and four- velocity vector of the matter 
distribution. Hereafter, the semi-colon denotes covariant differentiation. In 
general, £ is a function of time. 

The field equations in normal gauge for Lyra's manifold, as obtained by Sen 2 
are 

13 3 

R ij ~ 2 9i J R + ™ \9ii^ k = -SttGT^, (4) 

where <p is a time-like displacement field vector defined by 

& = (0,0,0,/3(t)), (5) 

and the other symbols have their usual meaning as in Riemannian geometry. 
The energy momentum tensor T lJ is not conserved in Lyra's geometry. Here we 
want to mention the fact that the ansdtz choosing the coordinate system with 
matter requir the vector field happens to be in the required form exactly in the 
matter comoving coordinates. The essential difference between the cosmological 
theories based on Lyra geometry and the Riemannian geometry lies in the fact 
that the constant vector displacement field p arises naturally from the concept 
of gauge in Lyra geometry whereas the cosmological constant A was introduced 
in adhoc fashion in the usual treatment. 

The field equations with the equations J2J) and (J5J take the form 

B B 2 B' 2 3 
2 B + Bt-A^ + ^ 2 = -™ (6) 
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The energy conservation equation is 



XP + -0$ + 







a- 





(10) 



where x — 87rG. Here and in what follows, a prime and a dot indicate partial 
differentiation with respect to x and t respectively. We note that the coefficient 
of bulk viscosity £ does not appear explicitly in the field equations above. For 
the specification of £, we assume that the fluid obeys a barotropic equation of 
state 

P = 1P, (11) 

where 7(0 < 7 < 1) is a constant. 

3 Solutions of the field equations and discussion 

On integrating the equation l(Tjl. we obtain 

A= T J < 12 > 

where I is an arbitrary function of x. Using equation (|12|l . equations © and 
JSJ) can be written as 

l-f9=0. (13) 

Since A and B are separable functions of x, so, ^- is a function of x. Conse- 
quently, equation l|13fl gives after integration 

B = lS(t), (14) 

where S(t) is an arbitrary function of t. Using the equation (|14fl . (|12fl becomes 

A = l jS. (15) 

The metric then takes the form 

ds 2 = dt 2 - S 2 (t) [dX 2 + e 2X (dy 2 + dz 2 )} , (16) 
where X = In /. Equations I© and give 

1 „ S S 3 _ n . _ . 

™=si- 2 s-si-V 3 ' (17) 

(18 » 
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Using the equation (|11[> and eliminating p(t) from equations $T7$ and (JTSJ, we 

have 

2§ + (l + 3 7 )|J - (1 + 3 7 )i + j(l - 7)/? 2 = 3 X £§- ( 19 ) 

In most of the investigations involving bulk viscosity, the coefficient of bulk 
viscosity is assumed to be a simple power function of the energy density 

m = toP n , (20) 

where £o and n are constants. If n = 1, equation (|20|l may correspond to a 
radiative fluid. However, more realistic models are based on n lying in the 
regime < n < i. 

Using equation (|2"U|> . (|T§|| can be written as 

2§ + (1 + 3 7 )y - (1 + 3 7 )^ + |(1 - 7 )^ 2 = 3 X £op"§. (21) 

In order to obtain the explicit dependence of S(t) on time, we assume the 
deceleration parameter to be constant, i.e., 



SS H + H 



S 2 \ # 2 



b (constant), (22) 



where H — ^ is the Hubble parameter. It is worthwhile to note that the 
solutions arrived at and the problem scanned cannot be properly handled by 
taking the DP to be a variable one. By considering DP to be a function of time 
i, we observe that q is rapidly increasing which is inconsistent and hence we 
should restrict ourselves to CDP. The equation (|22|l can be rewritten as 

■5+^=0. (23) 
Integration of equation 123[) gives us the exact solution 

S(t) = [a(t-t )]~^, when b ^ -1 (24) 
= m^™ 2 ', when b = -1 

where a, mi and mi are constants of integration and the constant to allows us 
the freedom of choosing the initial time. 
On using equation l|2"3)l in lt2"T)l . we obtain 

[3(1 + 7 ) - 2(1 + b)\ H 2 -(l + 3 7 )^ + |(1 - 7)/? 2 = ^ p n H. (25) 

Equation i|25[) with equation i|18|) leads to 

[3(1 + 7 ) -2(1 + ft)] fl 2 - (l + 3 7 )^ + ^(l- 7 )/3 2 



(26) 



G 



In what follows, we will solve the equations for two particular values of n, viz., 
n = and n = 1. 

3.1 Solutions for £ = £ 

In this case, we assume n = in equation (|2(J[1 . Equations l|2U|) and (|26l) become 
£ = £0 = constant and 

[3(1 + 7 ) - 2(1 + b)} H 2 - (1 + 3 7 )^ + 1(1 ~ 7 )/? 2 = SxtoH. (27) 

Case (i): b ^ — 1. For singular models, 5(0) = and hence, equation (|24|l leads 
to 

S = atT^T). (28) 

In the figure 1, we have shown how the scale factor evolves with time for two 
choices of b, b — 1 and b = 2. 




t 



Figure 1: The behaviour of scale factor for MODEL I (b ^ —1) with time for 
b = 1 (solid line) and 6 = 2 (dashed line) for a = 1. 



From equation (J2EJl, we have 



H = (TTbji' (29) 
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The evolution of the Hubble parameter with time is shown in figure 2 for 6=1 
and 6 = 2 for a = 1. 




Figure 2: The behaviour of Hubble parameter for MODEL I (b ^ —1) with time 
for b = 1 (solid line) and 6 = 2 (dashed line). 

Using the two equations (|28|l and (|29|l in the equation (|27(1 . we obtain 



/3 2 = 



3(i- 7 )(i + 6) 2 ; 2 



3x£o(l + &)t-(37- 26+l) + 



(3 7 + l)(l + 6) 2 

2b 

a 2 i < 1 + b > 



Using equation ijHUI) in l(T%)l . we get 
1 



XP 



(l- 7 )(l + 6) 2 t 2 



[2(2-6)-3 X &(l + 6)t] 



(30) 



(31) 



(1 - 7)a 2 £( 1 + b > 

The geometry of the universe, in this case, is described by the line-element 

ds 2 = dt 2 - aHT^t) [dX 2 + e 2X (dy 2 + dz 2 )} . (32) 
In equation l|3U|) . we observe that 1 = when t = t c , where the critical time t c 



satisfies 



[(37-2&+l) + 3 7 £o(l + &)*c] = 



(3 7 + l)(l + 6) 2 




-25 j 

6 ^0 

time 



Figure 3: The plot of energy density for MODEL I (b ^ — 1) with time for 6=1 
(solid line) and b — 2 (dashed line). (Here, we have taken x — 1> 7 = 0.4 and 
6 = 0.5) 



It is also seen that (3 2 > provided t > t c and (3 2 < provided t < t c . 

From equation l|31|). it is also observed that p > provided t > T c , where T c 
is the critical time given by 

T c [2(2 - b) - 3x^(1 + 6)T C ] = [ 2 1 ■ 

The figure 3 shows the evolution of the energy density with time for two values 
of b, viz., b — 1 and b = 2. We see from the figure that when t > T c , the energy 
density is positive, but when t < T c , the energy density becomes negative and 
then evolves to a constant negative value. 

Physical behaviour of the model: 

In the case of a non-flat model, when 6^—1, the Ricci scalar becomes 

R = — -p. ; - ■ (33) 

In equation (|33|l . we see that when t — > 0; (i) i? — > oo if 6 = 0, and (ii) 
R — * oo if b > 1. The equation (|33|l also suggests that when £ — > oo, i? — > if 
6 > 0. The expansion and shear scalars are given by 

9 = JYThjl' ^ 
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respectively. 

The model has a singularity at t = 0. As t — » oo, the expansion ceases. In 
this model, 5 = 0, which confirms the isotropic nature of the spacetime which 
we have obtained in 132(1. 

Case (ii): 6 = — 1. In this case, we obtain from the equation 1221) 

H = 0^H = H o = constant, (35) 
and scale factor is given by 

S = m ie m2t . 

The evolution of the scale factor with time in this case for mi = 1 and m-i = 2,3 
is shown in figure 4. 
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Figure 4: The behaviour of the scale factor for MODEL I (b = —1) with time 
for TO2 = 2 (solid line) and — 3 (dashed line). 



On using equation l|35|l . equations l|27|l and (|18fl reduce to 



/3 2 = 



3(1-7) 



3x^0-3(1+7)^ 



(1 + 37) -w' 



XP = 



(1-7) 



6H - 3x£oH 2 e " 



(36) 



(37) 



The geometry of the Universe, in this case, is described by the line element 

ds 2 = dt 2 - a 2 e 2Hot [dX 2 + e 2X (dy 2 + dz 2 )} . (38) 
From equation (|36|l . we see that 
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Figure 5: The behaviour of the energy density for MODEL I (6 = — 1) with time 
for m2 = 2 (solid line) and TO2 = 3 (dashed line). (Here, we have taken x = 1, 
mi = 1,7 = 0.4, Co = 0.4 and H = 0.5) 



1. /3 2 > provided Co < for all t > T, 

2. /3 2 > provided Co > ^f 1 for all t < T, 

3. (3 2 < provided Co > ^ for all t > T, 

4. /3 2 < provided Co < ^ for all i < T, 

5. /3 2 = provided f = T, 
where 

(1 + 37) 



T = 



2m 2 



In 



m 2 {3(l + 7 )ff 2 - 3 X £o#o} 



From equation (|37|l . we observe that 



2H 



(i) p>0, if Co < for all t > T 3 



2iH 

(ii) p>0, if Co > for all < t < T x 



X 
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where, 



T\ = — In 

2m 2 



In the figure 5, we have shown the evolution of the energy density with time 
for two values of to 2 , viz., to 2 = 2 and m 2 = 3. We see from the figure that 
the energy density increases with time and then approaches a constant positive 
value. 

Physical behaviour of the model: 

The Ricci scalar R is given by 

R = 2H o-^ETf (39) 
a e 

We clearly observe from equation l|39|) that (i) when t — > 0, R — > (2Hq — -i?), 
and (ii) when t — > 00, R — > 2iJp . The expansion and shear scalars are 

= 3Ho, a = 0. (40) 

This model represents an uniform expansion as can be seen from equation (|4L)|) . 
The flow of the fluid is geodetic as the acceleration vector /j = (0, 0, 0, 0). 

3.2 Solutions for £ = £ P 

In this case, we choose n = 1 and hence, equations (121 III and 1)26(1 now reduce to 

£ = Cop, (41) 

[3(1 + 7 ) - 2(1 + b)} (1 + 3 7 )i + ~(1 - 7)/? 2 = 

3£o# (3tf 2 -jg- \f) . (42) 

respectively. 

Equation l|42|l can be written as 

/?2 = 4 9g ff 3 - (3 7 - 2b + l)g 2 + (1 + 3 7 - 9£og)^ 2 
3 (1 - 7 + 3£ -ff) 

Case (i): b ± -1. 

Here we are considering a singular model. Hence, using equations i|28|) and l|29|l 
in equations l|43|l and (|18fl . we obtain 

/( 

P 2 = 



3a 2 (l + &) 2 t 2 [(l- 7 )(l + ^ + 3£ ] 
[9a 2 £ - 9f (l + bftT^) - (3 7 - 26+ 1)(1 + b)a 2 t + (1 + 3 7 )(1 + 6) 3 t^FF] (44) 

YO = 2 [(2-6)a 2 -2(l + b) 2 ^] 

X ^ a 2 (l + 6)i [(l_ 7 )(l + 6) t + 3e ] ■ 1 j 

Therefore, in this case, from equation l|44|) . we see that 
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time 



Figure 6: The behaviour of the energy density for MODEL II (6^—1) with 
time for b = 0.5 (solid line) and b — l(dashed line). (Here, we have taken X = 1, 
a = 1, 7 = 0.4 and £ = 0.4) 

1. f3 2 = if t = t c , 

2. f3 2 > if t = t c , 

3. /3 2 < if t < t c , 

where, the critical time t c satisfies the relation 

(1 + b) 2 tf™ [9^0 - (1 + 3 7 )(1 + b)t c ] = 9a 2 £ - (37 - 26 + 1)(1 + b)a 2 t c . 

It can also be seen from equation H45[) that p > if b < 2 with t > 0. The 
behaviour of the energy density with time is shown in the figure 6 for b = 0.5 
and 6=1. We see from the figure that the energy density approaches a small 
positive value after initially becoming negative for a short period of time. The 
physical properties of this model are similar to those discussed for model I. 
Case (ii): b = —1. 

As in the previous model, we now have 

H = Hq = constant. (46) 
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Equations (|43[l and l|18l) become 



4 [9£ m?ff 3 -3(l + 7)ff, 



(l + 3 7 -9^ -ffo)e- 2m2t ] 



XP = 



(l- 7 + 3Co#oK 
In equation l|47J) . we observe that 



[1 - 7 + 3^ ffo] 
[(37 - 1)6 



-2m 2 t 



37m?F 2 ] 



(47) 
(48) 



1.9998 



1.9996 



a 

T3 



1.9994 



w 1.9992 



1.999 



1.9988 




Figure 7: The behaviour of the energy density for MODEL II (b = — 1) with 
time for m 2 = 2 (solid line) and m 2 — 3 (dashed line). (Here, we have taken 
mi = l, 7 = 0.4, Co = -0.6, and H = 0.5) 

1. (3 2 >0 if Co < for all t < T 2 , 

2. /3 2 > if Co > for all t > T 2 , 

3. 1 < if Co < for all t > T 2 , 

4. /3 2 < if Co > for a11 * < T2 ' 

5. /3 2 = if Co ^ for all t = T 2 , 
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where 



T 9 



1 



2m 2 



-In 



(1 + 3 7 - 9£o#o) 



3m?flg(l+7-%Ho) 



From the equation Q48)). we see that 

(i) p > 0, provided £ > 
(ii) p > 0, provided £ < 



(7-1) 
3# 

(7-1) 
3# 



for all t > T 3 , 
for all i < T 3 , 



where 



T 3 



1 



2m 5 



-In 



37-1 
3 7 ^i^o 2 



The behaviour of the energy density with time is shown in figure 7. We see 
that the energy density increases with time and approaches a constant positive 
value. The physical properties of this model is similar to those of model I. 



4 Discussion and Conclusion 

In this paper, we have investigated LRS Bianchi type I models with a bulk vis- 
cous fluid and obtained exact solutions for a constant deceleration parameter. 
We have assumed the coefficient of bulk viscosity to be of the form = £o/? n ; 
where p is the energy density and n is the power index. The behaviour of the 
displacement field (3 and the energy density have been examined for values of 
n = and n = 1 for both (i) power-law and (ii) exponential expansion of a 
non-flat universe. The model discussed here is isotropic and homogeneous and 
in view of the assumption of isotropy, the shear viscosity cannot exist. The 
effect of bulk viscosity is to introduce a change in the perfect fluid model. 
Recently there is an upsurge of interest in scalar fields in general relativity and 
alternative theories of gravitation in the context of inflationary cosmology [SB] - 
Therefore the study of cosmological models in Lyra geometry may be relevant 
for inflationary models. Further the space dependence of the displacement field 
(3 is important for inhomogeneous models for the early stage of the evolution 
of the universe. Besides, the implication of Lyra's geometry for astrophysical 
interesting bodies is still an open question. The problem of equations of motion 
and gravitational radiation need investigation. Finally in spite of very good pos- 
sibility for Lyra's geometry to provide a theoretical foundation for Relativistic 
Gravitation, Astrophysics and Cosmology, the experimental point is yet to be 
undertaken. But still the theory needs a fair trial. 
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